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I review the use of the 2PI effective action in nonequilibrium quantum field theory. 
The approach enables one to find approximation schemes which circumvent long- 
standing problems of non-thermal or secular (unbounded) late-time evolutions en- 
countered in standard loop or 1/iV expansions of the IPI effective action. It is 
shown that late-time thermalization can be described from a numerical solution of 
the three-loo p 2EJ. effective a ction for a scalar <^ 4 -theory in 1+1 dimensions (with 
Jiirgen Cox, hep-ph/0006160). Quantitative results far from equilibrium beyond 
the weak coupling expansion can be obtained from the l/N expansion of the 2PI 
effective action at next-to-lea ding order fNLO) . calculated for a scalar O(N) sym- 
metric quantum field theory ( hcp-p h /01 0531 1 ). Exte nd inpr recent calculatio ns ir 



classical fiel d theory by Aarts et al. (hep-ph/0007357 ) and by Blagoev et al. (hep 
ph/0106195| ) to N > 1 we show that the NLO approximation converges to exact 
(MC) results already for moderate values of N (with Gert Aarts, hep-ph/0107129). 
I comment on characteristic time scales in scalar quantum field theory and the ap- 
plicability of classical field theory for sufficiently high initial occupation numbers. 



1 Introduction 

Nonequilibrium dynamics of quantum fields is much less well understood than 
its thermal equilibrium limit. Important progress has been achieved for sys- 
tems close to equilibrium, including (non)linear response techniques, gradient 
expansions or effective descriptions based on a separation of scales in the weak 
coupling limit EJ. Current and upcoming heavy ion collision experiments pro- 
vide an important motivation to find controlled approximation schemes which 
yield a quantitative description of far-from-equilibrium phenomena. Other 
applications include (p)reheating at the end of inflation in the early universe, 
the physics of baryogenesis or measurements of the dynamics of Bose-Einstein 
condensation. 

To study nonequilibrium phenomena no other dynamics than the known 
Hamiltonian time-evolution of quantum fields is necessary. In particular, the 
observed macroscopic, effectively dissipative behavior or thermalization has 
to arise from the underlying reversible quantum dynamics. Addressing these 
questions one has to find approximation schemes which respect all symmetries 
of the underlying theory, such as time-reflection symmetry. Far from equilib- 
rium the dynamics can involve very different length scales and a priori there 
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is typically no effective description based on a clear separation of scales which 
is valid at all times. Similarly, the relevant renormalized interactions between 
modes can become large at intermediate times and expansions based on weak 
coupling arguments are limited. 

Here we report on recent progress to find suitable approximation schemes 
for far-from-equilibrium dynamics based on^he two-particle irreducible (2PJ) 
generating functional for Green's functions 0. The approach does not employ 
the loop expansion of the 2PI effective action relevant at weak couplings aa, 
and extends first calculations of late-time Jjhermalization for quantum fields 
solving the three- loop 2PI approximation uu. To obtain a small nonpertur- 
bative expansion parameter for far-from-equilibrium dynamics we consider a 
systematic 1/N expansion of the 2PI effective actionH. Since the trun- 
cation error may be controlled by higher powers of 1 /N this expansion is not 
restricted to weak couplings or situations close to equilibrium. The next-to- 
leading order (NLO) approximation has been solved numjarically for a scalar 
O(N) symmetric quantum field theory in 1+1 dimensions Bill. The 1/7V expan- 
sion of the 2PI effective action in the NLO approximation is equivalent ±o the 
"BVA" approximation^] discussed earlier for quantum oscillators in Ref. 0. Re- 
cently, this approximation has been studied for a one-component classical field 
theory in Ref. □ and compared with exact (MC) results by numerically inte^ 
grating the classical evolution equations and sampling over initial conditionsEU 
Extending this analysis to the TV-component classical O(N) model one can 
showQ that the 1/N expansion at NLO yields systematically, by increasing N, 
quantitatively precise results already for moderate values of N. 

Controlled computational methods for the approximative solution of non- 
equilibrium dynamics are limited so far. Leading order large- N or mean field 
type approximations neglect scattering and are known to fail to-describe ther- 
malization E9. Standard (finite- loop) perturbative descriptions L!J or 1/N ex- 
pansions of the generating functional for one-particle irreducible {IP I) Green's 
functions beyond leading order can be secular (unbounded) and break dowja 
at late times It is known that the calculation of transport coefficients Ilia 
involves already an infinite series of diagrams not included at any finite order 
in the 1/N expansion of the 1PI effective action. In contrast, the relevant 
diagrams can be efficiently captured in finite-order loopt3 or 1/N expansion El 
at NLO and NNLO of the 2PI effective action. The three-loop order has been 
also frequently used as a^taxljpft^jpiut for kinetic descriptions and comprises 
the Boltzmann equation NoliMED. 

a The name "Bare-Vertex-iApproximation" refers to an ansatz on the level of a specific 
Schwinger-rDyson equation Q. Note that the four- vertex gets renormalized at next-to-leading 
order 2PIB. 
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Figure 1: Examples for the time dependence of the equal-time propagator G(t,t;p) with 
Fourier modes |p| = 0,3,5 from a first 2PI three-loop solution 4 (all in initial-time mass 
units). The evolution is shown for three very different initial conditions with the same 
average energy density. After an effective damping of rapid oscillations at early times the 
modes are still not close to equilibrium. It is followed by a smooth "drifting" of modes 
and a subsequent late-time approach to thermal equilibrium. The late-time behavior of the 
correlator modes is insensitive to the initial conditions and their value uniquely determined 
by the conserved energy density. This qualitative behavior has been verified beyond the 
loop/weak coupling expansion using the (2PI) 1/N expansion at next-to-leading order 6 . 



2 Thermal fixed point 

We consider a scalar quantum field theory with action, 

S = J d d+1 x (~p x a<p a d x oip a - ^d x (p a d x (p a - ^m 2 ip a tp a - (ip a <p a f ) (1) 

where (p a (x), x = (x , x), is a real, scalar field with a = 1, . . . , N components. 
For N = 4 this corresponds to the linear sigma model for the light scalar 
and pseudoscalar pion degrees of freedom (a, ir) respecting global chiral 0(4) 
symmetry, or to the Higgs sector of the electroweak standard model in absence 
of gauge field interactions. To save computational time the numerical results 
will be presented for 1+1 dimensions. We stress that in the low dimensional 
theory many realistic questions cannot be addressed, in particular since there 
is no spontaneous symmetry breaking. The dynamics may help to understand 
the symmetric regime relevant for sufficiently high energy densities. 

Before describing the nonequilibrium dynamics of thiSpmodel in more de- 
tail below we give here an overview. As pointed out in Ref.Q, for very different 
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noncquilibrium initial conditions a universal asymptotic late-time behavior of 
quantum fields can be observed. This behavior is demonstrated in Fig. |l] for 
the equal-time propagator G(t, t; p) in the three-loop 2PI approximation as a 
function of time t for three Fourier modes, each starting from three very differ- 
ent initial conditions with the same average energy density. For the solid line 
the initial conditions are close to a mean field thermal solution, the initial mode 
distribution for the dashed and the dashed-dotted lines deviate more and more 
from thermal equilibrium!]. It is striking to observe that propagator modes with 
very different initial values but with the same momentum p approach the same 
large-time value^j. The late-time result approaches the equilibrium correlation 
modes of the corresponding thermal field theory 

The observed real-time dynamics of correlator modes bears some simi- 
larities with a Wilsonian renormalization group flow of couplings towards an 
attractive fixed point. For the latter the large-distance or low-energy behavior 
near the fixed point is insensitive to details at short distances or high energies. 
For the real-time evolution we observe that the large-time behavior is insensi- 
tive to details at early times. The role of the fixed point is played here by a 
time-translation invariant solution for the correlators obeying the standard pe- 
riodicity (KMS) condition for thermal equilibrium. Time translation invariant 
solutions play an important role for the late-time dynamics of nonequilibrium 
field theory. This fact is known aspypllJbr the leading order (LO) large-iV 
or mean field type approximations E3E3o which neglect scattering: At LO 
the correlator modes-apoppach time-translation invariant solutions at asymp- 
totically large times EJEatil. Those leading order "fixed points'Ea are distinct 
from the thermal one. In particular, the LO fixed points depend explicitly 
on the initial particle number distribution, the simple reason beeing that the 
LO approximation exhibits an additional conserved quantity (particle number) 
which is not present in the full theory. The LO fixed points become unstable 
once scattering is taken into account and the system approaches the thermal 
fixed (cf. the detailed LO/NLO comparison in Ref. til). We emphasize that the 
described qualitative aspects are the same if we take into account scattering in 
the (2P/) three-loop or 1/N expansion at NLO, though quantitative aspects 
are different. We now go through the different characteristic time regimes in 
more detail using the 1/N expansion of the 2PI effective action at NLO, which 
provides for sufficiently large N a small nonperturbative expansion parameter. 

k We emphasize that thermalization or independence of the initial conditions cannot be ob- 
served in a strict sense because of time-reversal invariance (note that thermal equilibrium 
is invariant under time translation and bears no information about initial conditions). Our 
results demonstrate that thermal equilibrium can be approached closely with time, indistin- 
guishable for practical matters. It cannot reach it on a fundamental level without some kind 
of coarse graining which is a matter of principle and not a question of approximation. 



4 



3 Controlled nonperturbative approach 

All correlation functions of the quantum theory can be obtained from the effec- 
tive action T[(f>, G], here the 2PI generating functional for Green's functions, 
which is parametrized by the field G ab (x, y) representing the expectation value 
of the time ordered composite T(p a (x)<fif.(y) and the macroscopic field </> a (x) 
given by the expectation value of ip a (x)a We will concentrate in the following 
on the symmetric regime where it is sufficient to consider T[<fi = 0, G] = r[GL 
The 2PI generating functional for Green's functions can be parametrized asH 

T[G] = ^TrlnG" 1 + ^Tr G^G + T 2 [G] + const, (2) 

where Gq 1 = £(□ + to 2 ) denotes the free inverse propagator. WritingXa [G] = 
r^°[G] + r^ LO [G] + . . . the LO and NLO contributions are given byffl 

r^°[G] = ^d d+1 xG aa (x,x)G bb (x,x), (3) 

lf L °[G] =\J c d d+1 x ln[B(G)](x,x). (4) 

Here C denotes the Schwinger-Keldysh contour along the real time axis@ and 

B(x,y;G) = S^ +1 {x ~ y) + i-^ G ab (x,y)G ab {x,y). (5) 



- 1 . Li 
Ibyi 

ST[G] 



The nonlocal four-point vertex at NLO is given by ^ B 1 . £ n In absence of 
external sources the evolution equation for G is determined by 



5G ab (x,y) 



0. (6) 



In the following we solve Eq. (§), using Eqs. (§)-(!) 

without further approx- 
imations numerically in 1+ 1 dimensions. For a detailed description of the 
approach, the numerical implementation and results see Ref.Q. We present the 
results using the decomposition identity 

G(x, y) = F(x, y) - (i/2)p(x, y) sign c (x° - y°) (7) 

where F is the S5rmmetric or statistical two-point function and p denotes the 
spectral functional. In thermal equilibrium F and p would be related by the 
fluctuation-dissipation theorem, however, far from equilibrium both two-point 
functions are linearly independent. 
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\p(x,0;p=0)\ (NLO) 
\F(x°,0;p=0)\ (NLO) 




Figure 2: Logarithmic plot 6 of the spectral and statistical two-point modes, \p(x°, 0;p = 0)| 
and \F(x°, 0; p = 0)| (in units of Mi). The correlation modes oscillations quickly approach 
an exponentially damped behavior with a corresponding characteristic time scale. 



4 Early-time exponential damping 

We study the time evolution for two classes of initial condition scenarios away 
from equilibrium. The first one corresponds to a "quench" where initially at 
high temperature one considers the relaxation processes following an instant 
"cooling" described by a sudden drop in the effective mass. The second scenario 
is characterized by initially Gaussian distributed modes in a narrow momen- 
tum range around ipts , which is reminiscent of colliding wave packets moving 
with opposite and equal momentum c3. A similar nonthermal and radially 
symmetric distribution of highly populated modes may also be encountered in 
a "color .glass condensate" at saturated gluon density with typical momentum 
scale PtsEj- 

We start with a quench from an initial high temperature particle number 
distribution n (p) = l/(exp[ v /p 2 + Af 2 /T ] - 1) with T = 2M h Af 2 = 2M 2 
in units of the renormalized mass Mi at initial time t = 0, and a weak effective 
coupling A/67V = 0.083 M : 2 for N = 10. From Fig. | one observes that the 
statistical and spectral two-point functions, F and p, oscillate with initial 
frequency en/27r = 0.17Mt. The shown correlation functions quickly approach 
an exponentially damped behavior with characteristic rate ^ damp - ) — 0.016A/t. 
After the damping time scale r' damp - ) ~ l/7Q damp ' ) correlations with initial 
times are effectively suppressed and asymptotically F and p(t, 0;p) — > + . 
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Figure 3: Nonequilibrium damping rates extracted from F(t, 0;p = 0) as a function of 1/7V 
for strong initial coupling A/Aij 2 = 30. The quantum results are shown with full symbols. 
The damping rates are reduced compared to the classical field theory limit 7 (open symbols). 
For the classical theory round symbols represent results from the next-to-leading order l/N 
approximation and square ones denote the exact MC results. One observes a rapid conver- 
gence of the l/N expansion at NLO to the exact result already for moderate values of N. 



In Fig. H we show the damping rate as a function of l/N, this time for a 
strong initial coupling A/A/j 2 = 30. The initial particle number na(p) is chosen 
to represent a Gaussian distribution peaked around pts = 2.5A/i with a thermal 

[HI 



background of temperature To = 4Mi 



The full symbols show the NLO 



results for the quantum theory. One observes that the damping time ~ 1 /^ dam P 
scales proportional to N for sufficiently large N. It becomes apparent that LO 
(N — > oo ) or mean field type approximations fail to describe damping correctly 
and can be valid only for times t <C r( damp ). To demonstrate the power of the 
{IP I) l/N expansion we repeat the NLO calculation for the corresponding 
classical statistical field theory where comparison with exact (MC) results 0, 
including all orders of l/N, is possible. The systematic convergence of the 
NLO and the Monte Carlo result is apparent from Fig. ||. 

5 Non-exponential/power law drifting at intermediate times 

After the characteristic early-time scale r' damp - ) ~ l/^/ dam P the system is typ- 
ically still far away from equilibrium for a relatively long time. For a large 
variety of initial conditions we find a subsequent parametrically slow, non- 
exponential "drifting" of modes, as exemplified in Fig. [j]. To become more 
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Figure 4: Double logarithmic plot 6 for the equal-time mode F(x°,x°;p) for an initial particle 
number distribution peaked around p = pts ■ The dotted straight line corresponds to a power 
law behavior ~ (x° Minit) 1,75 ■ At later times the behavior changes to an exponential 
approach to thermal equilibrium with rate ^ thorm -ydamp 

quantitative we show in Fig. |i| the statistical two-point function on a double 
logarithmic plot for a Gaussian initial particle number distribution similar ta 
the one above and with small effective coupling X/6N = O.lA/j 2 for N = 1 0. Q 
Time-averaged over the oscillation period the evolution of the equal-time mode 
is well approximated by a power law behavior for times t < 30/A/tnit- The 
observed scaling behavior is reminiscent of a cascade in fully developed tur- 
bulence. However, note that the initial condition is homogeneous in space 
and all the energy is initially concentrated in a small momentum range. The 
energy has then to be distributed from the densely populated modes at high 
momentum in particular to the low momentum modes in order to reach ther- 
mal equilibrium. The presence of power law behavior at intermediate times 
can be observed for a variety of initial conditions. For higher n-point functions 
one can also observe approximate scaling, however, the detailed time evolution 
before the late-time thermalization can be rather complex and depends on the 
initial conditions. See Ref.Q for more details. 

6 Late-time exponential thermalization 

We emphasize that irrespective of the details of the initial conditions we find! 
quite different rates for damping and for the late-time exponential approach to 
thermal equilibrium. For the example of Fig. the characteristic time scales 
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Figure 5: Time dependence of the inverse slope T(t,p), as defined in the text 7,6 . When 
quantum thermal equilibrium with a Bose-Einstein distributed particle number n(f . e p ) is 
approached, all modes have to get equal T(f,p) = T cq . In contrast, for classical thermal 
equilibrium the inverse slope is momentum dependent with T(pi ow ) > T(phigh)- 



^(damp)^(thorm) ^ 0(io).| Though the exponential damping at early times is 
crucial for an effective loss of details of the initial conditions — a prerequisite 
for the approach to equilibrium — it does not determine the time scale for 
thermalization. It is possible to relate the damping rate to the "width" of the 
Wigncr transformed spectral function p in the same way as discussed in Ref. 
however, the thermalization rate cannot. This fact may be particularly pro- 
nounced in 1+ 1 dimensions where on-shell two- to- two scattering is constrained 
by the energy conservation relation such that it does not change the particle 
numbers for the involved momentum modes. Processes with nontrivial momen- 
tum exchange are necessary for thermalization and are taken into account at 
NLO by off-shell effects. Correspondingly the calculation of the three-partick 
peak from off-shell decay of the nonequilibrium spectral function in Ref. □ 
clearly shows a nonvanishing contribution. On-shell particle number changing 
processes appear at NNLO and may change quantitative aspects. The on-shell 
corrections beyond NLO are known to be small in the statistical field theory 
limit where exact results including all orders in 1 /N show very goad agreement 
with the NLO approximation already for moderate values of N M 



c For the employed spatially homogeneous initial conditions we do not observe a power-law 
"tail" for the late-time evolution. This behavior is insensitive to changes of large volumes. 
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The nontrivial fact that quantum thermal equilibrium is approached be- 
comes more pronounced if one compares with classical til thermalization. The 
classical field approximation is expected to become a reliable description for 
the quantum theory if the number of field quanta in each field mode is suffi- 
ciently high. Accordingly, we observe that increasing the initial particle num- 
ber density leads to a convergence of quantum and classical time evolutionLj. 
However, since classical and quantum thermal equilibrium are distinct the re- 
spective evolutions have to deviate at sufficiently late times, irrespective of the 
initial particle number density. Differences in the particle number distribu- 
tion can be conveniently discussed using the inverse slope parameter T(t,p) = 
—n(t, e p )[n(t, e p ) + \]{dn/ de)^ 1 for a given time-evolving particle number dis- 
tribution n{t, e p ) and dispersion relation e p (i)El Following Ref. tl we define the 
effective particle number as n(t, e p ) + i = [F(t, t';p) d t d t >F(t, t'\p)Y^ 2 \t=t' and 
mode energy by e p (t) = [d t d t 'F(t,t']p)/F(t,t';p)} 1 ^ 2 \t = t', which coincide with 
the usual free-field definition for A — > 0. For a Bose- Einstein distributed parti- 
cle number the parameter T(t,p) corresponds to the (momentum independent) 
temperature T(t,p) = T eq . In the classical limit the inverse slope T(t,p) as 
defined above remains momentum dependent. In Fig. [5] we plot the function 
T(t,p) for pi ow ~ and Phigh — 2pts- Initially one observes a very different 
behavior of T(t,p) for the low and high momentum modes, indicating that 
the system is far from equilibrium. Note that classical and quantum evolution 
agree very well for sufficiently high initial particle number densitytD. However, 
at later times the quantum evolution approaches quantum thermal equilibrium 
with a constant inverse slope T = 4.7MiEl In contrast, in the classical limit the 
slope parameter remains momentum dependent and the system relaxes towards 
classical thermal equilibrium u. 

7 Conclusions 

The 2PI effective action provides a powerful framework to find practicable 
approximations for nonequilibrium dynamics. Combined with a (2PJ) l/N 
expansion beyond leading order one obtains a controlled nonperturbative de- 
scription of far-from-equilibrium dynamics at early times as well as late-time 
thermalization. This includes the important capability of the I/TV expansion 
to describe nontrivial scaling properties near second order phase transitions at 
NLO and beyxmd. Most pressing is the extension to 3+ 1 dimensions along the 
lines of Ref. 13. This allows one to quantitatively address in a quantum field 
theory the formation of disoriented chiral condensates or fluctuations near crit- 
ical points in the context of heavy-ion collisions, or of (p)reheating at the end 
of inflation in the early universe. 
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